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I. INTRODUCTION 



There are two classical laws for the molecular liquids known more than a century. Long 
time they were considered as mere curious facts restricted to the simple van der Waals 
equation of state (vdW EoS). First of them is the law of the rectilinear diameter (LRD) 
]E 0] . It states that the diameter of the coexistence curve in terms of density-temperature 
is the straight line: 

ni + n a T r — T 

n « = J ^f = 1 + A ^rT • A>0 « 

where rij, i = l,g are the densities of the liquid and the gas phases correspondingly, n c is 
the critical density, T is the temperature. Further we put the Boltzmann constant to unit: 
fee — 1- Another simple linear relation is the Batschinski law || for the vdW EoS: 

P = -^--an\ (2) 
1 - no 

where P - is the pressure. It states that the line determined by the condition Z = 1 , where 
Z = P/(n T) is the compressibility factor is the straight line: 

n T 

- + ^ = 1- (3) 

Here n* = 1/6, T* = a/b is the Boyle temperature in the van der Waals approximation and 
a, b are the parameters of the vdW EoS (see e.g. Q). In general case of the spherically 
symmetrical potential T* is determined in accordance with H as following: 

T ^w) = a (4) 

where 



2vr / <$> attr {r)r 2 dr (5) 



and §attr{ r ) is the attractive part of the full potential $(r), a is the effective diameter of the 



particle so that b = ^ a 3 . The definition for the density parameter n t will be given below. 
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The connection between two linear relations |T] and (||) has attracted attention after work 
of Herschbach&Coll. |6|, 0, where the line Z = 1 was named by the Zeno-line and it was 



shown that it is indeed almost straight for the normal fluids though the deviations are 
noticeable for the accurate data. The straightness of the Zeno-line implies the constraint 
on the contact value of the radial distribution function || . Moreover, there are correlations 
between these linear elements and the locus of the critical point which were discovered in 
series of works of Apfelbaum and Vorob'ev P-[TT|]. The authors put forward the creative 
idea about the tangency of the Zeno-line to the liquid-vapor binodal extrapolated to the 
nonphysical region T — > 0. This allowed to connect the locus of the CP with the parameters 
of the Zeno-line: 

T B n B 

which are T B and n B . They determine the intersection points for (||) with the corresponding 
axes. According to former treatments || 0, 0, |10| , T B is the Boyle temperature, i.e. 

B 2 (T B ) = 0, 

and n B is determined by the relation: 

(7) 



T B dB 2 
UB ~B 3 (T B ) dT 



T=T B 



where B n is the virial coefficient of n-th order 



The relation ([!]) fails at the vicinity of the critical point (CP) where the singular terms 



appear Hl2| , |13| (see also recent works HUHlql). Nonlinear deviations from the linearities 



(H) and ^| in low temperature region, e.g. for water can be associated with the presence 



of the anisotropic interactions like H-bonds [17]. Despite the deviations from the exactly 



linear behavior the relations p] and (||) are much more general and can be considered as the 
nontrivial extension of the principle of corresponding states JF]. 

In |L8] it was shown that both p] and @ as well as the correlations between these linear 



elements and the locus of the critical point can be considered as the consequence of the 
global mapping between the liquid-vapor part of the phase diagram and that of the Lattice 
Gas. The latter is described by the Hamiltonian: 

H = ~J^2 qiqj - h ^ qi . (8) 

(ij) * 

Here qi is the site filling number and qi = 0, 1 whether the site is empty or occupied cor- 
respondingly. The quantity J is the energy of the site-site interaction of the nearest sites i 



and j, h is the field conjugated to the filling variable We denote by t the temperature 
variable corresponding to the Hamiltonian The density parameter is the probability of 
occupation of the lattice site x = ( ) . 

If the LRD is assumed then such mapping is given by: 

n=n *ITIt' T = T * TT~zt ' (9) 

where z, n* and T* are some parameters, which are connected with the coordinates of the 
CP: 

" (10) 



" c = 2(lt7)' 

T C =T*^—. (12) 
* 1+z K J 

This transformation is uniquely determined by the correspondence between the characteristic 
linear elements on the phase diagrams of the fluid and the LG. It is assumed that the 
coordinates of the CP for the LG are normalized so that t c = 1 and x c = 1/2. In such a 
context the parameter z represents the class of the corresponding states ||19|| . Note that at 
,2—t-O and T* — > oo with zT» — > 1 from |9| we get n/n* — » x and T — > t. 

The application of the transformation to the calculation of the locus of the critical point 



of Lennard- Jones fluids is given in [19, 20]. The strong argument in favor of the choice 



of the linear element 13 with the parameters n^T* instead of those for the Zeno-line 



and the comparison of the values for various potentials are given in [|19|, |2T|. There |9| was 
used to map the binodal of the planar Ising model onto the binodal of the two-dimensional 
Lennard- Jones fluid. 

The aim of this paper is to discuss physical basis of the linearities ([l]) and @ on the 
liquid- vapor part of the phase diagram of the fluids. We will follow the results of Jl8, |l9j. We 



expand some arguments of (jnj, which concern the corrected interpretation of the Batchinski 
law in a way consistent with the van der Waals approximation for the EoS. Also we derive 
the relation between the thermodynamic potentials for the continuum and the lattice models 
of fluids. 
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II. LIQUID- VAPOR BINODAL AS THE IMAGE OF THE BINODAL OF THE 
LATTICE MODEL 



It is easy to see that linear laws [5] and (^j) are fulfilled trivially in the case of the Lattice 
Gas (LG) model or equivalently the Ising model. Indeed, the rectilinear diameter law for the 
LG is fulfilled due to the symmetry of the Hamiltonian (|§|) with respect to the line x — 1/2. 

The analog of the Zeno-line for the LG can be defined too. In this case it is the line 
x = 1 where the "holes" are absent. This is consistent with the basic expression for the 
compressibility factor JIJ: 

Z = — = l-— r 3 ^r^g 2 (r;n,T)dr, (13) 



nT 3T J dr 

and the definition of the Zeno-line as the one where the correlations determined by the 
repulsive and the attractive parts of the potential compensate each other. It is clear that if 
x = 1 then the perfect configurational order takes place and the site-site correlation function 
of (H) vanishes: 

( ( ft Qj ) ) = ( ft Qj ) ~ ( ft ) ( Qj ) = • 

Thus the line x = 1 plays the role of the Zeno-line on the x — t phase diagram of the 
LG. Due to simple structure of the LG Hamiltonian (|J) and explicit symmetries there are 
degenerate elements of the phase diagram. The critical isochore x c = 1/2 coincides with the 
diameter. The Zeno-line is the tangent to the binodal which in this case expands up into 
the region t — > 0. Thus there is the degeneration of these in the case of the LG. Naturally, 
that both mentioned degenerations for the linear elements [I] and @ of the phase diagram 
disappear for the real fluids. The difference of the diameter and the isochore is nothing but 



the asymmetry of the binodal [22]. The difference between the Zeno-line and the tangent to 



the extrapolation of the binodal into the low temperature region T — > takes place only for 



the vdW EoS and influences directly the approach of |10], ^1] • The latter is based heavily of 
the constraint of the tangency to the extrapolation of the binodal. In fact there is no physical 
reasons to identify such tangent line with the Zeno-line. This is possible only for the vdW 



EoS. Using the generalized van der Waals approach of |24j any EoS can be approximated 
by the vdW EoS with the corresponding parameters. The definition of the tangent linear 
element ^] relies upon the van der Waals approximation for the given EoS and therefore does 
not coincide with the Zeno-line. The value of is determined by the condition analogous 



to ©: 

T* dB 2 
n *~B 3 (T*) dT 



(14) 

T=T* 



This relation follows from the constraint: 

d f Z(n(T),T)-l \ ( . 

df[ — W) — J {) 

which generalizes the condition Z = 1 and implies the linear change of the compressibility 
factor with the temperature T along the linear element (H). 

The proposed simple relation (Q) between the LG and the fluid may be useful for con- 
struction of the empirical EoS for the real substances. Recently, the nonlinear generalization 



of (Bp has been proposed in |25 



n=n *TT^' T=T *IT^' (16) 

with 7 as the fitting parameter. Unfortunately, the physical meaning of the parameter 7 
and its connection with the interaction potential was not discussed. Moreover from the 
basic thermodynamical reasonings the extensive parameters such as number of particles in 
the LG and in fluid should be proportional unless the fluctuational effects are taken into 
account. The lasts are the source of the fluctuational induced shift of the mean- field position 
of the critical point [EB . So the modification of the (^) in order to obtain the exact position 
of the critical point should be based on the inclusion of the fluctuation effects and their 
scaling properties. In J19 it was shown how the relations || augmented with some scaling 
considerations allow to obtain the critical points of the Lennard- Jones fluids basing on the 
properties of the potentials. The difference between Tb and T*, n# and n* is indeed essential 



especially in high dimensions d > 3 ||20|| . 

In the next Section we propose the generalization of the transformation (^) for the pro- 
cedure of the symmetrization of the binodal of the fluid which preserve the correspondence 
between the extensive thermodynamic quantities. 



III. SYMMETRIZATION OF THE BINODAL 

We believe that the key point which determines the form of the mapping is the LRD ([[]). 
The linear element |3] plays auxiliary role and defines the proper scales for the density and 
the temperature. 
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As was shown earlier in |21 the transformation (|9|) allows to map the binodal Tbi n (n) of 
the Lennard- Jones fluid onto the binodal £fem(x) of the corresponding lattice gas model. The 
latter has explicitly symmetric shape with respect to the critical isochore x c = 1/2: 



hin(x) = t bin (-x) , x = 1/2 — x, x < 1/2 (17) 

due to the particle-hole symmetry of the Hamiltonian §. Thus one can treat (0) as the 
procedure of the symmetrization of the phase diagram. Indeed, suppose that nd(T) is the 
dependence of the density diameter. Then the variable n/rid(T) is symmetrical over the 
binodal. With this © may be generalized as following: 

71 

x = ^n> r/r. = i - /(*) , (is) 

2n d {T) 

where the parametrization function / is chosen so that to map the binodal of the fluid 
Tbin{n) onto the binodal of the LG £&i n (x). It can be found from the common conditions 
of the thermodynamic equilibrium: P(n g ,T) = P(ni,T). The equality of the chemical 



potentials is fulfilled due to symmetry of the binodal of the LG (see also Section IV). 
The transformation (|9|) is the particular case for which: 

71 I 

n d (T) = ^f(t) 1 /(f) 



2 J w ' J w 1 + zt 

In the simplest linear approximation for the temperature behavior of the diameter of 
the form [l], obviously: 

n/n* = xf{t), T/% =1 -/(*). (19) 
The unknown parametrization function / is determined by the condition 

P(n h T)=P(n g ,T) . (20) 

To illustrate this procedure let us consider the classical vdW EoS (0). Using the symmetry of 
the LG representation |19] with respect to x c = 1/2 it is convenient to represent the densities 
of the coexisting phases as n g = x f(t) and rii = (1 — x) fit), < x < 1/2. Substituting these 
relations into (|20|) where the pressure P is given by (0) we get simple algebraic equation for 
the value /: 

l-2/ + / 2 + / 3 (x-l)x = 0, (21) 



as a function of x. In accordance with ( |T9D this provides the symmetrization of the binodal of 
the vdW EoS in terms of the LG variable x The result is shown in Fig. ID. The corresponding 



coordinates of the critical point in accordance with [19] are: 



n c /n, = -/(!) = 0.352, T c /T. = 1 - /(l) « 0.296 , P c /(n*T*) « 0.037. 



(22) 



The difference between the exact values n c /n t = 1/3, T c /T* = 8/27, P c / (n^T^) = 1/27 is 
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FIG. 1. The symmetrization of the vdW binodal using the parametrization (|l9|) corresponding to 
the linear approximation for the diameter. 



caused by the deviation of the diameter for the vdW EoS form the linear behavior. As we 
see the differences are rather small. Thus simple parametrization (|i~9|) can be applied to any 
EoS where the deviation from the LRD can be neglected. In this case the symmetrization 
of the binodal described by the parametrization (|I8~D or [19| along with the EoS for the LG 
could be useful in processing the data of the simulations for the coexistence curve (CC) of 
the Lennard- Jones fluids, where simple approximate formula: 



n, 



+ A\t\±B \t\P , 



T 



T-T c 



(23) 



is used p7| . Also the proposed approach allows to avoid ambiguity in extrapolation of the 
binodal into the region T — > fLOf . 

In order to visualize how the splitting of the degenerate elements mentioned above (see 
p. H) occurs it is expedient to consider the application of the transformation to the classical 



S 



EoS for the LG the Curie- Weiss molecular field approximation (see e.g. f2"8|| ). This equation 
of state has the formQ: 

h(m, t) = t ArcTanh (m) — m , m = 2x — 1 . (24) 

The isotherms for this model EoS are shown in (a). The linear character of the transfor- 
mation (|9]) with respect to the order parameters x and n allows to connect the pressure for 
the asymmetrical liquid Pu q (T,fi) with that of the LG as follows: 

P Hq (T, f i)=P LG (t(T),h( f t,T)). (25) 

The dependence h(fi,T) will be obtained in Section [TV| (see ([32|) below). The pressure for 
the LG Plg and in particular for EoS |4] can be defined by the standard method (see e.g. 



29fl). Taking into account the conjugation of the thermodynamic variables we can write: 



rn _|_ 

PLG(t,h) = f{t,m) + —^-h, (26) 



where f(t,m) is the thermodynamic potential for the variables (t,m). It is constructed 
easily using EoS (p4j) so that: 

F) f 

(27) 



nyt, m) 



dm t 

Substituting ( 2D ) into ( |2^| ) we are able to construct the isotherms in n — T plane. The 



result of construction of the isotherms and the binodal of the "fluid " in coordinates "pressure- 
density" basing on the Curie- Weiss EoS fl24]) for the LG is shown in ^|. 

In the following Section we discuss the relation between the thermodynamic functions of 
the LG and its continuum analog in detail. 

IV. THE RELATION BETWEEN THE THERMODYNAMIC POTENTIALS OF 
THE FLUID AND THE LATTICE GAS 



Basing on the analysis of |2T], [25| one can expect that the transformation ([|) gives the 
relation between phase diagram of the LG and that of the real fluid at least as the "zeroth 
order approximation". The idea that the difference between irregularity configuration for 



1 Here we neglect the trivial difference between the m-field conjugated to variable m and the x-field conju- 
gated to x because it is irrelevant for our consideration. 
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0.4 0.6 

(b) 

FIG. 2. The isotherms obtained using @ applied to Curie- Weiss approximation for LG with 
z = 1/2. The isotherms of (^) are shown in Fig. (a). The corresponding isotherms for the 
isomorphic "fluid" model are in Fig. (b). The red curve is the binodal, the critical isotherm 
T c = 1/3 is marked by the blue line. The quantities are dimensionless (see (|22|) ). 



continuum fluids and regularity of configurations of lattice models is unimportant for con- 
sideration of the order-disorder transitions in the fructuational region is due to K.S. Pitzer 
(see P0|). But beyond the fructuational region the shape of the holes in real or continuum 
liquid and that in the lattice gas causes the main difference between the configurations of 
these systems. From this point of view the global character of the transformation |9| shows 
that the particle-hole simplified picture still can be useful. Though, it is not the particle 
density which reflects such symmetry. Rather the combination of the density of the particles 
and the density of the holes is the symmetrical variable (see p. |13|). This rehabilitates the 
hole theory for expanded liquids |29|, j3l] . Such caricature picture of the liquid state gives 
the possibility to relate the thermodynamic functions of these systems. 
Let 

<5(t,h,J\r)=ArQ(t,h), and J(T,fi,V) = P(T,fi)V , 

are the thermodynamic potentials of the grand canonical ensembles for the LG and the fluid 
correspondingly. Here M is the number of sites in a lattice. First it is natural to state the 
following relation M = n* V between the extensive variables of these ensembles. The results 



of [32] allow to interpret 1/n* as the volume per particle in the ideal crystal state at T — > 0. 
Such a state defines the lattice which may serve as the basis for the determination of the 
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corresponding Lattice Gas model. 
Using the standard definitions: 



1 ^1 



1 d<5 



x 



M dh 



(28) 



along with @ we get the following relation between the potentials: 

J(/i, T,V) = <5{ h{n, T),t(T),N) P(/i, T) = n* ( T),t(T) 



(29) 



From (P) the relation between the density of the fluid and the density of the LG can be 
written as following: 



n(//, T)/n* = x(h(ji, T),t{T)) ( 1 - T/T t ) . 
Taking into account trivial relation: 



(30) 



d 




dh 


d 


d ji 


T 


d n 


t dh t 



from (f28"l) ,f29| and ( p0|) we get the following relation: 

h( f j,,T) = (l-T/T,)(fi-f M) (T)) , 



(31) 



or, in the inverse form: 

At-A* (T) = /i (1 + zt) . (32) 

We remind that h = below CP is the coexistence line for the LG and is mapped onto 
the saturation curve of the continuum fluid. Therefore imq(T) coincides with the chemical 
potential fi s (T) along the saturation curve below the critical point T < T c . To determine 
(io(T) in the supercritical region T > T c we note that the line 

1 - T/T* 



nGuo(T), T)/n* 



(33) 



is the image of the line of symmetry h = for the LG along which x = 1/2. Therefore Ho{T) 
defined by (|33|) can be considered as the Widom-Stillinger line of symmetry for the liquid 
|33| . Indeed, such line is defined as the locus of maximum of correlation length where the 
thermodynamic response function such as the isothermal compressibility Xt — - jt has 

TL Op rp 

maxima T max (P) j34j- Obviously, in case of the LG it is represented by the line h = 0, or 
equivalently by the critical isochore x — 1/2. 
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FIG. 3. Widom line (dashed) as the line of the temperature maxima of the isothermal compress- 
ibility xt(P,T) of the fluid for different P. It is given by (|33|) . The compressibility is calculated 
using the Curie- Weiss EoS (p3) with subsequent transformation (pi). 



The Widom-Stillinger line has attracted much attention recently in the studies on super- 
critical states of liquids |J5[ . The stated mapping between the LG and the fluid naturally 
explains the separation of the fluid region into the gas-like and the liquid-like states with the 
Widom line as the border. The gas-like region is the image of the region x < 1/2 where the 
number of holes is greater than the number of particles. For the liquid-like region x > 1/2 
the situation is inverse. 

The line ([J3|) is the image of the line h = 0. Therefore along this line the continuation 
of the subcritical behavior (divergence) of the isothermal compressibility \t takes place. 
Fig. |] shows the result of the calculations for the compressibility xt in accordance with the 
relation ^between the fluid and the LG. This is the locus of the compressibility maxima. 

From (P^| ) the relation between the entropy & of the LG and that of the fluid S can be 
derived: 



S(T,fi,V) 



dJ_ 



dt d<5 



1 + zt 



e(t, h,M) 



(34) 



dT dt 

The results described above illustrate the general analysis of [13|. The mapping ^ deter- 
mines the fluid as the asymmetric model for the symmetric LG. The conservation of the RDL 
is due to the direct relation between the the thermodynamic potentials EM Indeed, if linear 
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relation between P(n, T) and g{h(fi, T), t(/i, T)) with the analytic coefficients holds, then the 



quantity 



dP(ii,T) 



dh 



has the rectilinear diameter. This order parameter is the combination of 



the density and the entropy O, [L5 



In order to be physically meaningful, the transformation |] or its generalization should 
appear as the average of the transformation of the microscopic variables. We assume that 
there exist such transformation which leads to the symmetrization of the binodal. The 
similarity between the liquid-vapor part of the diagram for the fluid and that of the LG and 
also the separation between the liquid-like and gas -like states (see e.g. f3lj ) allows to state 
that for the case of fluids there exists the "asymmetrical" microscopic observable A such 
that: 



[A) l + (A)=0. 



(35) 



Here {...) lg stands for the averages on the corresponding coexisting liquid and gaseous 
states with (T, ni(T)) and (T,n g (T)). In terms of |37j the equilibrium average (A) lg 7^ is 
the density-like variable, which takes different values in the coexisting phases. The specific 
form of the observable depends on the Hamiltonian of the fluid. Such observable exists in 
case of penetrable sphere models where the relation between the thermodynamic potentials 
P and g is linear [fl3|| . In fact the right hand side in ^5] can be any analytic function of T 
(including the neighborhood of the CP). In such case it is possible to redefine A so that ( |35"1) 
is fulfilled. Then the following equation 



[A) 







(36) 



determines the continuation fio{T) of the diameter into the supercritical region T > T c . Once 
such observable is determined the corresponding thermodynamic potential can be defined 
so that 



(A) 



dh A 



In general 5 depends on the field variables P, fi, T nonlinearly. Therefore, in accordance 



with [[L3|, [37| the density diameter shows both 1 — a and 2(3 anomalies. Since the singularity 
of the diameter is of fluctuational nature the transformation (|9]) should be considered as the 
mean-field approximation for the transformation between microscopic fields (see |16| ) . 
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V. THE CUBIC SHAPE OF THE BINODAL 



The fact of the global cubic shape of the CC for the molecular liquids is the long standing 
issue and was well known to van der Waals due to studies of Verschaffelt [[38| (see also the 
review [33]). It is also well established fact for a wide variety of the Lennard- Jones fluids 



with the short ranged interactions. The results of computer simulations are well described 
by the Guggenheim-like expression |23] (see 0). 

The transformation (||) states that the shape of the CC for the molecular fluids is de- 
termined by that for the LG [2D]. As is known from the computer simulations [[HJ the 



crossover to the classical behavior which is characterized by the parabolic shape of the bin- 
odal is absent for the LG with the nearest neighbors interaction. On the basis of the global 
character of the transformation (|9]) one can assume that the same is true for the molecular 
fluids with the short-ranged potentials. In other words, the binodal of the molecular fluid 



can be described approximately by the global cubic dependence similar to ( |23|) in a broad 
temperature interval. In |21| this statement was demonstrated for 2D and 3D Lennard- Jones 
fluids. In particular the binodal of the 2D Lennard- Jones fluid was obtained as the image of 
the binodal of the 2D Ising model given by the Onsager exact solution. In such a case very 
flat shape of the binodal dome of 2D fluid is due to the exponent (3 = 1/8 and the global 
nature of the power-like dependence of the binodal of the 2D LG: 

. = 1/2 ± /(*)>/• , /W-i-see^j- 

Note that f(t) is the analytic function of the LG temperature variable t. 

The classical Landau theory of the phase transitions gives the general picture which is 
characterized by the classical exponents for the critical asymptotics of the thermodynamic 
quantities ||. In particular the dome of the binodal is given by the the quadratic curve. 
But it should be stressed that the Landau theory breaks the continual character of change 
of the thermodynamic state in passing through the CP. Indeed, there is the discontinuity 
of the specific heat Cy and therefore there is no unique critical state but still there are two 
coexisting phases. This discontinuity is connected with the assumption of the existence of 
the only strongly fluctuating quantity. For the liquid-vapor critical point such quantity is 
the density. Indeed, if the entropy is used as the order parameter then the finite jump of 
the specific heat allows to distinguish between the ordered and the disordered phase at the 
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(a) LG (b) liquid 

FIG. 4. The phase diagram of the LG (the Ising model) and its mapping to (n, T) plane of liquid 
states with the help of @. The points mark the positions of the inflection points (see the text). 
The red point is the CP. Two families with h > and h < of iso-h lines t(h, m) (a) and their 
images in n — T plane (b) are shown. The binodal is the limiting curve of these families. 

CP. Besides, the parabolic shape of the binodal is the consequence of the usage of analytic 
EoS with applying the Maxwell construction for the pair of conjugated variables (e.g. the 
pressure and the volume) in the subcritical region. This puts the constraint of analyticity 
of the binodal in terms of other pair of the variables, e.g. the temperature and the entropy. 
Note that in the vicinity of the critical point t(x) should be the even function of the order 
parameter (see (|T7|) ). One can expect that weak divergence of the specific heat admits the 
analyticity of the function t(x). Then t(x) = 1 — box 2n + . . . , 6 > 0, where n > is 
the integer. Obviously, this occurs in d = 2 and d > 4, where the specific heat has weak 
logarithmic divergence or finite jump with n = 4 and n = 1 correspondingly. For d = 3 the 
specific heat diverges more strongly. Therefore one could expect that the analyticity of t(x) 
breaks in this case. 

To clarify where the possibility for the global cubic shape of the binodal appears let us 
consider the t — m phase diagram of the LG (Ising model). We use the Ising-like order 
parameter m = 2x — 1 for the convenience and consider the family of the iso-h curves 
t(h,m) for the LG. Notably, all these curves has the inflection points m (h). Fig. |] shows 
the situation for the Curie- Weiss EoS fl24|). Obviously, in the vicinity of the inflection point 
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m (h) determined by the condition: 

d 2 t(h, m) 



m=mo 



dm 2 

the function t(h,m) has the cubic form (see |j): 

t{h, m) = t (h) + h(h) ( m - m (h) ) + t -^- (m - m {h)) 3 + . . . , (37) 

o 

where the coefficients ti(h) and ts(h) are the odd functions of h. Naturally, the binodal as 
the line of the phase equilibrium is the curve which correspond to h = 0. From this point 
of view the binodal consists of two parts. They are the limiting curves of the families of 
cubics (|37|) with h > and h < correspondingly (see Fig. f|). Therefore the binodal as the 
limiting curve for these families should have the cubic form: 



£ (o) 

t(0 ± 0, m) = 1 - ^- m 2 \m\ + . . . . (38) 



or, equivalent ly 



where 



m = B Q \r\ 1/s + . . . , (39) 

B ° = (<?) 1/3 ' (40) 

In the mean-field approximation t 3 (h) — > oo in the CP. But the analysis of the simplest 
Curie- Weiss approximation shows that it is reasonable to use the extrapolated value tf^ 
(see Fig. [5]). From (flff ) we obtain the value of the amplitude for the binodal opening 
Bo 1.7. Surprisingly, this coincides with the result of the corresponding value obtained 
in the computer simulations for the 3D Ising model Under the transformation ^ the 

iso-h lines map onto the corresponding ones in n — T plane (see Fig. ^). Obviously, the local 
cubic character similar to fl3"T| ) near the inflection points takes place. Approaching the CP 
we obtain the relation between the amplitudes of the LG and the continuum fluid: 

= (l + Z y/ 3 B . (41) 



In accordance with the results [20 1 for the Lennard- Jones fluids, where the attractive part 
of the interaction has the "— 1/r 6 " behavior, z = 1/2 in d = 3. Then, from (40) and (|41 ) we 
get: 

B ( LJ) ^ 1.95, (42) 
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FIG. 5. The behavior of the coefficient b{h) in the Curie- Weiss approximation. The extrapolation 
is shown by the dashed line. 



which is good agreement with the value obtained in computer simulations Bq 
for potentials of "6 — n " type W2 . 



1.92 2.0 



This allows to treat the value (3 = 1/3 as the corresponding exponent which conforms 
with the continuity of the critical state for d = 3 for the systems isomorphic to the LG. Other 
critical indices, except the small ones a and rj which in the mean-field approximation are 
zeroth, can be obtained via the standard thermodynamic stability and scaling reasonings, 
so that v = 2/3, 7 = 4/3, 5 = 5 etc. 

VI. CONCLUSIONS 



In this paper the explicit relations between the basic thermodynamic functions of the 
LG and the continuum fluid is derived within the global isomorphism approach proposed 
in 



[L8| . These relations allow to obtain the information about liquid state directly from the 
EoS for the lattice model. It is quite remarkable since the lattice models do not contain the 



translational degrees of freedom. Nevertheless, as it follows from the results [21, 25] this 
mapping gives rather good description for the LJ fluids. In particular the stated relations 
can be used to connect the computer simulations for the lattice models with those for the 
continuum fluids (see [0]). Within such an approach the splitting between critical isochore 
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and the diameter, the Zeno-line and the tangent to the binodal at T — > is naturally 
described. The conservation of linear character of the diameter allows to obtain simple 
form of the transformation between the phase diagram of the LG and the continuum fluid. 
Obviously, the value of the critical density is the most sensitive to such approximation. But 
for the computer simulations the deviations from the law of rectilinear diameters are small 
to observe them near the CP P7[ . This explains rather good agreement of the estimates for 



the locus of the CP based on the parametrization (^j for the Lennard- Jones fluids with the 
results of the computer simulations [|19|]. 

We believe that the approach proposed could be useful in studies of the supercritical 
behavior basing on the EoS of the lattice models especially for water [34], [43| where a lot 



of lattice models, including the models with H-bond are known. In view of the results of 
work it would be interesting to generalize the proposed approach in order to search the 



correspondence between the dynamic response functions of fluids and their lattice analogues. 
This poses the question about the choice of the adequate lattice model of the proper geometry 
and the type of the interaction which realizes the isomorphism with the corresponding fluid. 
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